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We discuss a method to solve models with long-range interactions in the microcanonical and 
canonical ensemble. The method closely follows the one introduced by R.S. Ellis, Physica D 133, 
106 (1999), which uses large deviation techniques. We show how it can be adapted to obtain the 
solution of a large class of simple models, which can show ensemble inequivalence. The model 
Hamiltonian can have both discrete (Ising, Potts) and continuous (HMF, Free Electron Laser) state 
C^ . variables. This latter extension gives access to the comparison with dynamics and to the study 

jy^ ' of non-equilibrium effects. We treat both infinite range and slowly decreasing interactions and, in 

particular, we present the solution of the a-Ising model in one-dimension with < a < 1. 
C^ , Keyvirords: Long-range interactions. Large deviation techniques. Mean-field limit. 
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I PACS numbers: 05.20.-y Classical statistical mechanics 
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Q ' A system with long-range interactions is characterized by an interparticle potential V{r) which decreases at large 

O I distances r slower than a power r~" with a < d, d being the dimension of the embedding space [l|. Classical 
examples are self-gravitating |2| and Coulomb J3 systems, vortices in two-dimensional fluid mechanics [j|, wave- 
particles interaction 0, |g and trapped charged particles [3 • The behaviour of such systems is interesting both from 
J>. , the dynamical point of view, because they display peculiar quasi-stationary states that are related to the underlying 
00 ' Vlasov equations |8|, and from the static point of view, because equilibrium statistical mechanics shows new types of 
VO ! phase transitions and cases of ensemble inequivalence Q- In this paper we will restrict ourselves to the second aspect. 
C^ . In long-range interacting systems, essentially all the particles contribute to the local field: the fiuctuations around 

the mean value are small because of the law of large numbers. This explains qualitatively why the mean-field scaling, 
which amounts to let the number of particles go to infinity at fixed volume [l^, llil, is usually extremely good. 
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f-^ ■ However, let us remind that for long-range interac ting systems, microcanonical and canonical ensembles are not 
necessarily equivalent in the mean-field limit |l3.ll3.ll^. Moreover, because of the non additivity of the energy, the 

C^ ' usual construction of the canonical ensemble cannot be applied. This is the reason why the microcanonical ensemble 
is considered by some authors [^ Il5l | as the only "physically motivated" one. Hence, it is extremely important to 

T ' develop rigorous techniques to solve non trivial physical models in the microcanonical ensemble. One finds in books 
■ O ] the solution for the perfect gas, but generalizations to interacting particle systems are difficult. 

^ ■ The goal of this paper is to advocate the use of large deviation techniques as a tool to explicitly derive microcanonical 

and canonical equilibrium solutions for a wide class of models. As a first step in this direction, we will discuss here a 
general solution method to treat in full detail mean-field models without short distance singularity. Large deviation 
techniques |l6lll7j are nowadays widely used. For example, Michel and Robert |18J . using these techniques, rigorously 
derived the statistical mechanics of the two-dimensional Euler equations. The method was later used by Ellis et 

^ • al. [l9J to study other two-dimensional geophysical fluids. The statistical mechanics of some models with discrete 

Ci I variables has been recently obtained using large deviation theory |3, |2l|, |22| . Besides presenting the solution method, 
we will show in this paper how it can be applied to models whose Hamiltonian depends on continuous state variables, 
like the so-called Hamiltonian Mean-Field (HMF) model |2J,|23; the interest being here to obtain solutions in the 
microcanonical ensemble of models which display a Hamiltonian dynamics, opening the possibility of studying also 
non-equilibrium effects. 

We will first briefly introduce in Section Q] the large deviation technique. In Section^ we will recall the different 
steps of the mathematical framework introduced in Ref. [l^ necessary for a systematic application to long-range 
interacting systems: we will use the infinite range Potts model as a simple example. In a first instance, we will 
then treat infinite range models with continuous variables. We will present in Section IlII Al the solution of the HMF 
model [23 in the microcanonical ensemble and, in the following Section Fill Bl we will consider the Colson-Bonifacio 



model of Free Electron Lasers (FEL) [23|, as an example of relevance for physical applications. In Section ITVI we 
will show that these techniques can be applied also to cases where the interaction is not infinite range but distance 
dependent; the solution of the so-called a-Ising model |23, in one dimension with < a < 1, will be presented in full 
detail, allowing us to discuss also the role of boundary conditions. Finally, Section will be devoted to conclusions 
and perspectives. 

I. LARGE DEVIATION THEORY 

We will present in this section the main ideas behind large deviation techni q ues , with an emphasis to applications. 
For a more rigorous mathematical treatment, we direct the reader to Refs. |lfiLll7|. 

A. The large deviation principle 

Let us consider the sample mean of N independent real random variables Xk with the same distribution and zero 
average 

1 ^ 

fc=i 

The law of large numbers states that Sn tends toward the average x = (Xk), namely 0, when N tends toward infinity. 
Moreover, if Xk has a finite variance, since all hypotheses of the central limit theorem are fulfilled, the probability 
distribution P{\/NSn & [x,x + dx]) converges towards a Gaussian, hence the fluctuations of Sn are of order 1/\/N. 
Typical questions of the large deviation theory are: What is the behavior of the tails of the distribution? What is the 
probability of a fluctuation of order one of Sn, i.e. what is the value of P{Sn G [x,x + dx])? 

Let us be more specific by discussing the usual example of the coin toss. We attribute to heads and tails of a coin 
the values X^ = +1 and Xk = —I, respectively. The sum Sn can take {N + 1) distinct x- values in the interval [—1, 1]. 
For such values, using simple combinatorial analysis, one easily derives the probability distribution 



P{Sn^x)^ /,, , ^^,,x /,, ^^,,^ , (2) 

I! 2^ 



{l+x)N\, fil-x)N\,^j, 



which, using the Stirling's formula, can be approximated in the large N limit as 

In P{Sn = x) N f^i^^ In (I + x) + -^^^ \n{l - x) ] = -NI{x) , (3) 



which defines the function I{x). More precisely, one can prove that for any interval ]a;i, X2[G [—1, 1], 

lim \nP{x g]xi, 3^2 [) = max I{x) . (4) 

AT— >CX3 A^ x£]xi, X2[ 

In the language of large deviation theory, one states that Sn fulfills a large deviation principle, characterized by 
the rate function I{x). If one interprets the coin toss experiment as a microscopic realization of a chain of N non- 
interacting Ising spins, it is straightforward to prove that, in the statistical mechanics vocabulary, I[x) corresponds 
to the negative of the Boltzmann entropy (divided by the Boltzmann constant) of a state characterized by a fraction 
X of up-spins. This is a first simple example of the large deviation principle, and the main purpose of this paper is to 
present other examples of its use for more complicated and physically relevant systems. 



B. Cramer's theorem 

Cramer's theorem |l7| allows one to derive the probability distribution P{Sn G [x,x + dx]) in the large A'^-limit, 
providing also a method to eompute the rate function I{x). The theorem is formulated for multi-dimensional and 
identically distributed random variables Xk G M'', d being the dimension of the space of the variables. We will 
formulate the theorem in an informal way, without emphasizing mathematical technicalities. 

Let us define the function ^(A) as 

*(A) = (e^-^) , (5) 

where A G M'', "•" is the usual scalar product and ( ) is the average over the common probability distribution of 
the variables Xk- Cramer's theorem states that, if ^(A) < oo, VA G R'', then the sample mean Sn satisfies the large 
deviation principle 

\iiP{Sn e[x,x + dx])^ -NI{x) , (6) 

with rate function I{x) {x £ M.'^) given by the Legendre-Fenchel's transform of ln$, 

I{x) = sup (A-2;-ln1'(A)) . (7) 

The formulation of the theorem is not restricted to discrete random variables: this will be important for our applica- 
tions. 

A heuristic proof of the theorem for the simplest case Xk € M goes as follows. The probability of obtaining Sn = x, 
with d/i the common probability distribution of each variable Xk , is given by 

N 

P{Sn e [x, x + dx]) = H d^i{Xk) 5{Sn - x) . (8) 

■^ fc=i 

This formula can also be interpreted as the volume of the phase-space {Xi, ■ ■ ■ ,Xjv) under the microcanonical con- 
straint that Sn = x. Using the Laplace transform of the Dirac's (5-function, one obtains 

P{Sn e [x,x + dx]) = ^ d\ e-^^- / l[dfi{Xk) e^S"-i^'= , (9) 

^* -^^ -^ k=l 

where F is a path on the complex A-plane going from —ioo to -|-ioo, which crosses the real axis at a positive value. 
Subsequent manipulations of this formula lead to 

P{Sn e[x,x + dx]) = ^J' dX e-^^- [(e^^)] "" = ±-f d\ e-^(^--i"<-'")) "^-"^ e'^^^^) , (10) 

where I{x) is given in formula Q with d = 1. In the last step, a large N saddle-point approximation has been 
performed. 

Most of the results contained in this paper will be obtained using Cramer's theorem, because the statistical variables 
of the models we will consider are identically distributed in space (mostly on a lattice) . In all cases the function In ^P 
is differentiable, which also fulfils the hypotheses of the Gartner-Ellis theorem [l7l |. 

II. A GENERAL METHOD 

In this Section, we will describe the use of the large deviation method to solve models with long-range interactions. 
As already mentioned, Michel and Robert [l3| successfully used large deviations techniques to rigorously prove the 
applicability of statistical mechanics to two-dimensional fluid mechanics, proposed earlier |23,|33. Ellis et al. [13 have 



developed and generalized this approach to solve two-dimensional geophysical systems with long-range interactions. 
Here, we will adopt Ellis et aZ.'s approach, emphasizing the different steps in the construction of thermodynamic 
functions. The method will be exemplified discussing in detail the three-state Potts model with infinite range inter- 
actions. This simple example has been recently used as a toy model to illustrate peculiar thermodynamic properties 
of long-range systems |31| . The diluted three-state Potts model with short-range interactions has also been studied 
in connection with "small" systems thermodynamics by Gross |32l |. 
The Hamiltonian of the three-state Potts model is 

J ^ 
Hn^^^Y.^s,.s,. (11) 

The 1/N prcfactor is introduced in order to keep energy extensive jsj. Each lattice site i is occupied by a spin 
variable Si, which assumes three different states a, b, or c. A pair of spins gives a ferromagnetic contribution —J 
(J > 0) to the total energy if they are in the same state, and no contribution otherwise. It is important to stress that 
the energy sum is extended over all pairs («, j): the interaction is infinite range. 
The solution method consists of three steps. 

Step 1 : Identifying global variables 

Let S^v be the phase-space of a TV-particles system with Hamiltonian 

Hn ■■ Sat ^ M (12) 

and o^jv G Sat be a specific microscopic configuration. The first step of the method consists in associating to every 
microscopic configuration lon, a global (coarse-grained) variable, 7(a;Ar). Then, a new Hamiltonian H^ can be defined 

Hn{u!n) ^ Hn{j{u!n)) + RN{i^N) ■ (13) 

If one can neglect i?Ar(wAr) with respect to Hn in the large A^-limit, then the Hamiltonian can be expressed only in 
terms of the global variables. When considering the above defined infinite range Potts model, the appropriate global 
variable is 

7 = (na^nb^Uc) , (14) 

where (na,nb,nc = 1 — Ua — rib) are the fractions of spins in the three different states a,b,c. In this case, the 
Hamiltonian expressed in terms of the global variable is 

HN^-^{nl+nl+nl) , (15) 

and coincides with the original Hamiltonian Hn even at finite N. In Section llVI we will discuss the a-Ising model, 
for which Rn does not vanish. 

The global variable 7 is of finite dimension in our example, but could be of infinite dimension in other cases. For 
instance, 7 could correspond to a local mass density in a gravitational system, or a coarse-grained vorticity density 
in 2D turbulence. 

Although this type of redefinition of the Hamiltonian might appear to be possible in all cases, this is not true. 
For instance, in the case of short-range interactions, even after defining a local density of a physical quantity (e.g. 
magnetization), when performing the large TV- limit, no general argument exists to neglect Rn- As a consequence, in 
such a case, the local density (e.g. the local magnetization) is not the appropriate macroscopic variable. However we 
will argue that such a procedure is, instead, viable in general for systems with long-range interactions on a lattice. 
Keeping the lattice spacing finite is important in order to regularize possible short distance singularities. 

Step 2 : Deriving an entropy functional for the global variables. 



Because the mean-field variables are not equiprobable, the number of microscopic configurations leading to a given 
value of 7 does depend on 7 itself. Then, one can define an entropy functional 5(7), 

s(7)= lim llnf}Ar(7) = -/(7) + lnAA , (16) 

N~>oo I\ 

where the Boltzmann constant has been set to unity, Qn (7) is the number of microscopic configurations corresponding 
to a given value of 7, I{-f) is the rate function and A/" = / d'yV,]\[{-f) is the total number of states. This is where large 
deviation theory applies; not only in proving that such an entropy functional exists, but also in providing a procedure 
to derive it explicitly. 

If, besides the dynamical variables which contribute to the global ones, the Hamiltonian depends also on a number 
of variables n„ which is small with respect to N, it is easy to prove that these additional variables contribute to the 
entropy for a negligible term, proportional to ny/N. 

For the infinite range Potts model, it is possible to derive the entropy functional 5(7) directly, using combinatorial 
arguments and the Stirling approximation in the large iV-limit. However, let us follow instead the procedure given 
by Cramer's theorem as explained in Section llBl Expression (I14|l for 7, which is identified with x in Section llBl can 
be rewritten as 



(nYI ^s.m, ^ Y. ^s.,b, ^ J2 '^•5-- 1 



\ i i i 

The local random variables with common probability distribution are here 

Xk^{5sk,a,5su,b,5s^,,c) ■ (18) 

Hence, the generating function ^P is given by 

^'iXaAbAc) = ^ E (e^"*--+^^^-.^+^^^--) (19) 

= \ {e^- + e^' + e^^) . (20) 

o 

The large deviation functional is 

7(7)= sup {XaTia + XbUb + XcTlc - In *(Aa , Afc, Ac)) . (21) 

This variational problem can be solved exactly, giving \i = Inn^ with i ~ a,b, c. Hence, 

^(7) = nalnua + Ublnnb + {I- na-nb)ln(l -Ua-rib) + \ii3 . (22) 

Thermodynamic entropy density is given by s{j) = —I{j) + In A/", where the normalization factor is A/" = 3 for the 
Potts model example, which recovers the result of the combinatorial approach. 

Step 3 : Microcanonical and canonical variational problems 

To obtain the entropy as a function of energy density e, i.e. to solve a model in the microcanonical ensemble, after 
performing steps 1 and 2 of the solution method, one has to find the solution of the following variational problem [I^l 

5(e) = sup (s(7) I i/(7) = e) , (23) 

7 



where 



Hi,) = hm ^ . (24) 
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FIG. 1: Caloric curve (inverse temperature vs. energy density) of the three states infinite range Potts model. The canonical 
solution is represented by a solid line. The microcanonical solution coincides with the canonical one for e < e' and is instead 
indicated by the dash-dotted line for e* < e < —J/6. The increasing part of the microcanonical dash-dotted line corresponds 
to a negative specific heat region. In the canonical ensemble, the model displays a first order phase transition at /?'. The 
two dotted regions bounded by the dashed line and by the microcanonical dash-dotted line have the same area (Maxwell's 
construction). 



This result is exact but it corresponds also to an "intuitive" mean-field solution. Finally, let us notice that the total 
entropy density is intensive: there is no difference in this respect with short-range interacting systems. When other 
conserved quantities exist, they must be taken into account when solving the variational problem in formula ^2'6\ . 
Specific examples will be discussed in sections (IIII All 



and (|IirB|l . 

For the infinite range Potts model (|II|I . the variational problem is 



S{e) = sup i-nalnua - mliinb - (1 - «„ - n6)ln(l - Ua - nt) -— {nl + nf + {I - Ua - nb)"^) = ej 



(25) 



This variational problem can be solved numerically. The microcanonical inverse temperature (3{e) = dS/de can then 
be derived: it is shown in Fig. Q] in the allowed energy range [—J/2, —J/6]. Ispolatov and Cohen |31| have obtained 
the same result by determining the density of states. A negative specific heat region appears in the energy range 
[-0.215 J, -J/6]. 

Let us now consider the canonical ensemble. The partition function, written in terms of the common probability 
distribution of the phase-space variables, is 



ZW,N) 



N 

= / YidfiiXk)^ 
•' k=i 



'PH, 



N 



(26) 



where ^l is the probability density oi Xk- This is not the usual partition function, but differs from it only for a constant 
factor, which counts the number of states. Let us remark that we have also used the letter [3 for the microcanonical 



inverse temperature. We will comment specifically when the microcanonical inverse temperature differs from the 
canonical one due to ensemble inequivalence. 

Approximating Hamiltonian i/yv hi the large TV-limit with the one expressed in terms of the global variable 7, one 
gets 

/3i/(7) 

For infinite range models, formula (|27|l is exact for all N because the rest Rn vanishes. Using formula H16() . one 
obtains 

Z(/3,7V) c /d7e-^(-^(7) + /3^(7)) . (28) 



d-fnNh)e~ 
Z{f],N) '^^°° ^ . . (27) 



Applying the saddle point method, the partition function is rewritten as 

Z{f3,N) ~ e-^^('^) , (29) 

where the "free energy" F{p) is obtained solving the variational problem 

F(/3)=inf(/37J(7)-s(7)) ■ (30) 

7 

Our "free energy" is the usual free energy multiplied by the inverse temperature. This helps because physical states 
will correspond to minima of such free energy also for negative inverse temperatures. 

In the case of the infinite range three-state Potts model, the canonical free energy can be explicitly derived solving 
the following variational problem 

F{/3) ^ inf (nalnna + nblnnb + nclnric — (nl + nl + nl) \na + Ub + ric ^ l) ■ (31) 

To obtain the caloric curve, one has to compute e = dF / dj3. Fig. ^ shows that at the canonical transition inverse 
temperature /3* ~ 2.75, corresponding to the energy e*/ J ~ —0.255, a first order phase transition appears, with an 
associated latent heat. The low energy "magnetized" phase becomes unstable, while the high energy "homogeneous" 
phase, which has the constant energy density, e/ J = —1/6, is stabilized. In Fig. ^ the two dotted regions have the 
same area, respecting Maxwell's construction. At the inverse transition temperature, there is also a jump in the global 
variables (na,nf),nc) which are the order parameters of the model. 

This extremely simple example shows already ensemble inequivalence. In the microcanonical ensemble, there is no 
phase transition and the specific heat becomes negative. On the other hand, in the canonical ensemble, there is a first 
order phase transition with a latent heat. The caloric curves do not coincide. We observe that in the energy range of 
ensemble inequivalence, microcanonical temperatures, {dS/de)~^, do not coincide with any canonical one. 

This is an example of the more general fact that entropy s{£) is not always the Legendre-Fenchel transform of the 
free energy F[(3). A general discussion of ensemble i neq uivalence, both at the thermodynamic level and at the level 
of equilibrium macrostates, is provided in Refs. |l9l l22l |. while a classification of phase transition and of ensemble 
inequivalence situations is reported in Ref. |23l |. 

This concludes the general presentation of the different steps of the method to derive the statistical mechanics of 
long-range interacting systems. 

III. EXAMPLES 

In this Section, we will discuss the application of the large deviation method to two examples which share the 
difficulty of computing the entropy for a phase-space with continuous variables. When presenting the method, we 



have discussed in parallel its application to a model with discrete variables, the infinite range three-state Potts model, 
which however could have been solved by direct states counting J3l| . This latter approach cannot be used when the 
variables are continuous. Obtaining microcanonical entropy often implies the solution of too complicated integrals, 
and indeed one docs not find many examples of such solutions in the literature. On the contrary, the large deviation 
method is not restricted to discrete variables and we will show that it can even be simpler to use in such a case. 

A. The Hamiltonian Mean Field model 



The Hamiltonian Mean Field (HMF) model |24L l25l| is defined by the following Hamiltonian 

^^ = Ey + ^E^°«(^^-^^) ' (32) 

where di € [0, 27r[ is the position (angle) of the i-th article on a circle and "Pi the corresponding conjugate variable. 
This system can be seen as representing particles moving on a unit circle interacting via an infinite range attractive 
(C < 0) or repulsive (C > 0) cosine potential or, alternatively, as classical XY-rotors with infinite range ferromagnetic 
(C < 0) or antiferromagnetic (C > 0) couplings. The renormalization factor N of the potential energy is kept not 
only for historical reasons, but also because it simplifies the derivation of the variational problems and makes the 
problem well defined. As we will see, this implies that the usual energy per particle and temperature are well defined 
in the N —f oo limit. In the literature, some authors have treated the case in which the energy is not extensive. This 
leads to different thermodynamic limit behaviors |3J, |43| . 

The canonical solution of this model has been derived using the Hubbard-Stratonovich transformati on l25l| . The 
microcanonical solution has been heuristically obtained, under the hypothesis of concave entropy in Ref. |35l | and in a 
different form in Ref. [33 ■ In this section, we will derive both solutions with no additional hypothesis. We will verify 
that the two ensembles give equivalent predictions. 

Step 1: Hamiltonian (|32|) can be rewritten as 



^iv=V^ + ^(Mj + M2) ^ (33) 



^-^2 2 

=1 



where the magnetization is 

M = M^ + iMy = ^ X! ^''^'"' ' ^''^^^ ^'^ ^ 1^1 ■ (34) 

k 

By a direct inspection of Hamiltonian (|33|) . one can identify the global quantities u = ^ J^i Ph ^^^x and My. Moreover, 
since v = j, "YliPi i^ ^ conserved quantity with respect to the dynamics defined by Hamiltonian H32|) . it will be included 
in the global variable. Hence, 

-1={u,v,M,,My) . (35) 

The Hamiltonian in terms of the global variable is 

H{-i) = \{u + CM^) . (36) 

Step 2: The vector of local variables is Xk = {pi, Pk, cos Ok, sm0k) ■ The generating function is 

*(A„,A„,Ae,A,) = LXuP^ + XvP + Xc cos e + Xs sin e\ (37^ 

- e-^"/^^"y^/o(yAfTA!) , (38) 



where Iq is the modified Bessel function of order 0. In the last expression, we have not reported the constant factor 
Jy, dp dO which is finite because the domain of integration E is bounded due to the finiteness of the energy. The sign ~ 
indicates that formula (|38(l is valid only at leading order; indeed, as S 7^ M x [0, 27r[, neither the Gaussian nor the 
Bessel functions are fully exact. 

The large deviation functional is then given, apart from trivial constants, by 



/(7) 



sup 



\uu + Xyv + XcAU + XsMy - In/o (y^fTAf j + ^ + ^ ln(-A„) 



(39) 



This variational problem can be solved for the "kinetic" subspace (A„, A„) separately from the "configurational" one 
(Ac, As). The entropy as a function of the global variable is then 



■5(7) = Skin{u,v) + Sconf{M) 



where 



Skin{u,v) 



1 



In (m — w") + const 



^/(M) = -sup[AA/-ln/o(A)] = -A7l/ + ln/o(A) 



(40) 

(41) 
(42) 



with A = yAj+Aj and A the solution of the variational problem in 1)42(1 . Let us remark that Cauchy-Schwarz 
inequality implies that u> v^ . 



Step 3: It is therefore possible to derive the microcanonical variational problem 



S{e,v) 



sup 

M,u 



-ln(u-z;2) + Sco„/(M) 



i.+ fM- 


= e,u = 


= const 



and the canonical one 



F(/3) = inf 

M.u.v 



(3 



{u + CM^) - - ln{u - v') - s,onf{M) 



(43) 



(44) 



We show here the solution of both variational problems in the case v = and C = — 1, for which a second order 
phase transition appears. Shifting v to non- vanishing values does not produce anything new, whereas the second order 
phase transition disappears for positive values of C. 

Solving the sup condition in Eq. (|42|l leads to the consistency equation M = /i(A)//o(A) = m(A), which determines 



the optimal value of A, A 



^{M). Using the energy constraint, Eq. H43|l can be rewritten as 



Sis) 



sup 

Mg[04[ 



s(Af,e) = -ln(2e + A/2) 



^conf 



(M) 



(45) 



In order to determine the global maximum of S, let us note that in Eq. (|45l) M takes values in the interval [0,1[. 

Moreover, an asymptotic expansion shows that s{M, e) ^ (1/2) ln(l — M), which diverges to —00 when M tends to 
1. Thus, a global maximum of the continuous function s with respect to M exists and is attained inside the interval 
[0, 1[. To determine this maximum one has first to solve the extremal condition 



M 



2e + M 



5 £ 

conj 



(Af) 



-(M) 



(46) 



The unique solution of this equation is Af = for £ > 1/4, while a non vanishing magnetization solution bifurcates 
from it at e = 1/4, originating the second order phase transition. One can indeed show that 



s{M,e) 



ln(2e) 



Te-'^'' 



oiKP 



(47) 
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FIG. 2: Entropy versus energy (a) and free energy versus inverse temperature (b) for the HMF model (I33I I with C = — 1 
V — 0. The dotted hues are traced at the phase transition point. 



and 



which clarifies the stability change at e = 1/4 of the M — solution. Formula (|47|) was also obtained in Ref. |2g in 
connection with the fluctuations of the magnetization. 

Let us remark that {dS)/{de) is nothing but the inverse of twice the kinetic energy, which is the usual microcanonical 
temperature. Moreover, condition l|4()|) coincides, as expected, with the consistency relation derived in the canonical 
ensemble [23 ■ In Fig. [3 we show the full dependence of the entropy on energy and the graph of the free energy 
versus the inverse temperature, obtained by solving numerically the consistency equation (|46|) in the low energy 
(temperature) range. In Fig. |3 we plot the caloric curve and the dependence of the order parameter on energy: the 
two ensembles give the same predictions because the entropy is concave. 

In Appendix A, we discuss a different way of obtaining the microcanonical solution, which treats the kinetic part 
of the energy in a more traditional way, like, for instance, in gravitational dynamics |37l |. The method showed in this 
section is of more general applicability. 



The Colson-Bonifacio model for the Free Electron Laser 



In the linear Free Electron Laser (FEL), a relativistic electron beam propagates through a spatially periodic magnetic 
field, interacting with the co-propagating electromagnetic wave; lasing occurs when the electrons bunch in a subluminar 
beat wave |27l |. Scaling away the time dependence of the phenomenon and introducing appropriate variables, it is 
possible to catch the essence of the asymptotic state by studying the classical Hamiltonian 



N 



H 



N 



E? 



N5A' 



N 

2 A ^ sin(6'j 



-^) 



(48) 



The Pi's represent the velocities relative to the center of mass of the N electrons and the conjugated variables 6i 
characterize their positions with respect to the co-propagating wave. The complex electromagnetic field variable, 
A = Ae^"^, defines the amplitude and the phase of the dominating mode (A and A* are conjugate variables). 5 is a 
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FIG. 3: Inverse temperature versus energy (a) and magnetization versus energy (b) for the HMF model 13311 with C = — 1 and 
« = 0. The dotted hnes are traced at the phase transition point. 

parameter which measures the average deviation from the resonance condition. In addition to the "energy" iJ, the 
total momentum P = ^^Pj + NA^ is also a conserved quantity. Most of the studies of this model have concentrated 
on the numerical solution of Hamiltonian (|48|l . starting from initial states virith a small field A and the electrons 
uniformly distributed with a small kinetic energy. Then, the growth of the field has been observed and its asymptotic 
value determined from the numerics. Our study below allows to find the asymptotic value of the field analytically. 

Step 1: Similarly to the HMF case, Hamiltonian H48|) can be rewritten as 



H 



N 



NH{j) = N (--SA^ + 2A (-A4 sin ip + My cos (p)) 



(49) 



where M^, My, u and v have been defined in Eq. H34() and following. Defining the phase of the mean field ip' as 

Mx + iMy — M exp («(^'), the global variable is 7 = (u, u, M, (^', A, (p). 

Step 2: As remarked in Step 2 of Section |nl the contribution to the entropy of the two field variables A, ip, is 
negligible (of order 1/N). Hence, the VP function reduces to the one of the HMF model, see formula (|38|) . Finally, one 
obtains the same contributions to the kinetic and configurational entropies, as shown in formulas (|41|l and (|42(l . 

Step 3: Defining the total momentum density as ct = P/N, the microcanonical variational problem to be solved is 



S{e,a,S) ~ sup 



■ In (U- 1)2) + Scon/ (M) 



e = - + 2AM sin {ip' - p) - 5A^, a = v + A^ 



(50) 



Using the constraints of the variational problem, one can express u and v as functions of the other variables, 
obtaining the following form of the entropy 



S{e, a, S) = sup 

A,tp,M,{p' 



In 



2 e 



4AM sin {ip' -ip) + 2{6 - a)A^ - A^ 



./(M) 



(51) 
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The extremization over the variables ip and </?' is straightforward, since by direct inspection of formula (|51ll . it is clear 
that the entropy is maximized when ip' — tp = —-n 12. Then 



S'(e, CT, S) = sup 



AM 



i.n 



(7^ 



2 ( e - — ) + 4AM + 2((5 - a) A' - A" 



.f{M) 



sup s{A, M) . (52) 



A,M 



The non zero a case can be reduced to the vanishing a problem using the identity 5(e, ct, 5) ~ S{e — cr^/2, 0, (5 — ct). 
From now on we will discuss only the zero momentum case, changing e + g^/2 — > e and (5 + ct — > (5. This has also a 
practical interest, because it is the experimentally relevant initial condition |38l |. 
The conditions for having a local stationary point are 



ds _ 2{6A~A^ + M) 



dA 2e + 2(5^2 + AAM - A^ 
ds _ 2A 

dM ^ 2e + 28 A^ + ^AM - A'^ 



, (53) 

m-i(Af) = , (54) 



where m"^ is defined in formula (|46|1 . It is clear that A/ = A = is a solution of conditions H53|l and H54(l : it 
exists only for positive e. We will limit ourselves to study its stability. It must be remarked that this is the typical 
initial condition studied experimentally in the FEL: it corresponds to having a beat wave with zero amplitude and 
the electrons uniformly distributed. The lasing phenomenon is revealed by an exponential growth of both A and the 
electron bunching parameter Af . 

The second order derivatives of the entropy s(A, Af), computed on this solution, are 

£(0.0)^1 |i(0.0)^-2, ^(0.0) ^i . (55) 

The two eigenvalues of the Hessian are the solutions of the equation 

5\ 25 1 



.^-.^-2+-j----.0 . (56) 

The stationary point is a maximum if the roots of this equation are both negative. This implies that their sum 
S = (—2 + S/e) is negative and their product V = —2S/e — l/e^ is positive. Recalling that we restrict to positive e 
values, the condition for the sum to be negative is e > S/2 and the one for the product to be positive is e > —l/{26) 
with S < 0. The second condition is more restrictive, hence the only region where the solution M ^ A = exists 
and is stable is e > —1/(2(5) with 5 < 0. When crossing the line e = —1/(25) {S < 0), a non zero bunching solution 
(A/ ^ 0) originates continuously from the zero bunching one, producing a second order phase transition. This analysis 
fully coincides with the one performed in the canonical ensemble in Ref. [39l |. 

The maximum entropy solution in the region complementary to the one where the zero bunching solution is stable 
can be obtained |38l | by solving numerically Eqs. H53(l and H54|l . This corresponds to having a non zero field intensity 
and bunching. 

We have not completed in this case the study of the global stability of the different solutions, but we think that, 
in view of the possibility to map this model exactly onto the HMF model (see Appendix B), no surprise is expected 
and that the study presented here should fully represent all physical solutions. 

IV. THE ISING MODEL WITH l/r° INTERACTIONS 

All the models that have been considered above are infinite range: interactions are independent of the distance and 
the energy can therefore be written exactly in terms of global variables. This is no more valid for several important and 
physically relevant cases. Let us mention in particular the 1/r interaction law for gravity and Coulomb systems and 
the logarithmic interaction for two-dimensional turbulence. A generalized 1/r" interaction has been also introduced 
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and the resulting phase transitions have been analyzed |40j. In all these cases the interaction is singular at short- 
range. Interesting toy models (which generalize the HMF model) without a short distance singularity have been 
recently proposed [411 [42, |43, 13 ■ I^^ these latter models, XY spins are put on a lattice and interact through a slowly 
decrea sing non integrable 1/r" law (a < d). For what concerns studies in the microcanonical ensemble, is has been 
shown |42| that the thermodynamic behavior of these models is independent of the a exponent, after one has adopted 
an appropriate renormalization of energy and temperature. Salazar et al. |4fil | have studied numerically the problem 
using microcanonical Monte-Carlo simulations. Besides confirming the scaling properties with a, they have also shown 
that the number of states is of order exp{N). The exact solution of these models in the canonical ensemble has been 
obtained by Campa et al. |43 and VoUmayr-Lee and Luijten [4J|. The scaling of the magnetization and of the energy 
curve with the a-exponent has been exhibited in full detail. 

In this section, we will present a microcanonical solution of the one dimensional a-Ising model, whose Hamiltonian 
is given by 

where J > and spins Si = ±1 sit on a one-dimensional lattice with unitary spacing. The N°'~^ prefactor is 
introduced in order to have an extensive energy. This model has been first introduced by Dyson J28l | and studied for 
the "integrable" case, a > 1, in the canonical ensemble without the N°'~^ prefactor. We will show that it is possible 
to obtain an exact microcanonical solution using large deviation theory, when < a < 1. This solution can be easily 
generalized to lattices of larger dimension. In a preliminary paper [46j , a microcanonical solution of this model was 
presented without rigorously proving the exactness of the mean- field limit. 

The study of this model will also give us the opportunity to emphasize the important role played by boundary 
conditions when the interactions arc long-range. We will indeed consider both free and periodic boundary conditions. 
In the latter case, \i — j\ is the minimal distance along the circle where one identifies the first and the (A^ + l)-th 
lattice point. 

In the solution wc will adopt the same scheme described in Section Hll 

Step 1: The Hamiltonian Hn cannot be rewritten exactly using a finite dimensional global variable. We overcome 
this difficulty by defining a coarse-grained function. Let us divide the lattice in K boxes, each with n = N/K sites, and 
let us introduce the average magnetization in each box ruk, k = 1 . . . K. In the limit N —^ oo, K —^ oo, K/N — > 0, the 
magnetization becomes a continuous function m{x), of the [0, 1] interval. After a long but straightforward calculation, 
described in Appendix C, we can show that it is possible to express Hj^ as a functional of m{x): 

Hn = NH[m{x)] + o{N) , (58) 

where 

The estimation is uniform on all configurations. 

Step 2: The probability to get a given magnetization m^. in the fc-tli box from all a priori equiprobable microscopic 
configurations obeys a local large deviation principle P{mk} oc CKp[ns{m.k)], with 

s(TOfe) = — In — — In — . (60) 

Since the microscopic random variables in the different boxes are independent and no global constraints has yet been 
imposed, the probability of the full global variable (mi, . . . , mj^) can be expressed in a factorized form as 



K K 

„ns{mi) 



P{m,,m,,...,mK) = H^^"^') =^ 11^"''^"''^ = ^^P 



^ K 



,NS[m{x)] ^ (-g;^-) 
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where S[m{x)] ~ L s{m{x)) dx is the entropy functional associated to the global variable m{x). Large deviation 
techniques rigorously justify these calculations |47l |. proving that entropy is proportional to N , also in the presence 
of long-range interactions. This result is independent of the specific model considered; it applies, for instance, also to 
the long-range XY spin model studied by Salazar et al. [i^l • 



Step 3: It is now straightforward to formulate the variational problem in the microcanonical ensemble 

5'(e) = sup {S[m{x)] \e = H[m{x)] ) . 

ni{x) 



(62) 



Let us remark that the optimization problem l|62() has to be solved in a functional space. In general, this has 
to be done numerically, taking into account boundary conditions. In this paper, we consider only free and periodic 
boundary conditions. In the former case, the only available solutions arc numerical. An example of a maximal entropy 
magnetization profile obtained for free boundary conditions is shown in Fig.0|for different values of a. In the following 
of this section, we will treat the periodic boundary conditions case, for which analytical result can be obtained. 
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FIG. 4: Equilibrium magnetization profile for the a-Ising model with free boundary conditions at an energy density e = 0.1 
for a — 0.2 (solid line), a = 0.5 (dotted line) and a = 0.8 (dashed line). 



In the periodic boundary case, the distance \x — y\ in the energy (|57|l is defined as the minimal one on the circle, 
obtained when the two boundaries of the interval [0, 1] are identified. Both entropy and free energy can be obtained 
in analytical form for homogeneous magnetization profiles. 

In Appendix D, we prove that for /3 < /3c = (1 — a)/( J2") there is a unique global maximum of S, corresponding 
to a constant zero magnetization profile. The variational problem H62|l . where S is defined in Eq. (|6()|l . leads to the 
equation 



tanh (771(2;)) ~ P J 



m{y) 







x- 



y[ 



dy 



(63) 



where /? is a Lagrange multiplier. For /3 > /3c, we restrict ourselves to solutions with constant magnetization profiles, 
i.e. m{x) ~ m. We prove in Appendix D that these solutions are locally stable, i.e. close non constant profiles have 
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ca 




FIG. 5: a) Equilibrium magnetization in the allowed energy range in the microcanonical ensemble for the Q-Ising model with 
a = 0.5; the negative branch is also reported with a dotted line, b) Inverse temperature versus energy in the microcanonical 
ensemble (solid line). The canonical ensemble result superposes to the microcanonical one in the interval [/3c, oo] and is 
represented by a dashed line for (3 G [0,/3c]. Pc is then the inverse critical temperature in the canonical ensemble. In the 
microcanonical ensemble, no phase transition is present. 



a smaller entropy. Moreover, these are the only solutions when a = 0, since the right-hand-side of Eq. H63() is then 
independent of x. For constant profiles, the relation between energy and magnetization is 



1 



(64) 



where we have used /„ dx\x — y|~" = 2"/(l — a) and Emax = l/(2/3c)- Hence, fixing the energy impUcs fixing the 
magnetization and. consequently, the Lagrange multiplier (3 in Eq. H63|) . Expressing the magnetization in terms of the 
energy in the entropy formula (|60|l allows to derive the caloric curve (3 = dS/de. The consistency equation H63|) has 
always a non vanishing magnetization solution in the whole energy range [0,eniax]: this is reported in Fig.[Sfa). The 
caloric curve is shown in Fig. El^b) with full line. The limit temperature (3c is attained at zero magnetization, which 
is a boundary point. 

In the canonical ensemble, one has to solve the variational problem H3U|) . This leads to exactly the same consistency 
equation H63|l . where the Lagrange multiplier is replaced by the inverse temperature (3. Solving this consistency 
equation on the full positive (3 axis, one finds a zero magnetization for (3 < (3c and non vanishing one for (3 > (3c- The 
caloric curve is obtained taking the derivative e = dF/d(3. The graph of this function is reported in Fig. Elb) and 
superposes to the microcanonical caloric curve from infinity down to (3c while it is represented by the dashed line for 
(3<(3c. 

It follows that in the region [0,/3c], the two ensembles are not equivalent. In this case, a single microcanonical 
state at emax corresponds to many canonical states with canonical inverse temperatures in the range [0,/3c[. Thus, 
in such a case, the canonical inverse temperature is not equal to the microcanonical one. In the microcanonical 
ensemble, the full high temperature region is absent and, therefore, no phase transition is present or, in other terms, 
the phase transition is at the boundary of the accessible energy values. The entropy is always concave, hence no 
inequivalence can be present in the allowed energy range, apart from the boundaries. This situation is called partial 
equivalence [l9ll20|. This ensemble inequivalence persists for all a values below one, and is removed only for a = 1 
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when Eniax —^ oo and /3c — > 0: the phase transition is not present in both ensembles and the system is always in its 
magnetized phase. 

The main drawback of this analysis is the difhculty to obtain analytical solutions of Eq. H63(l for non constant 
magnetization profiles, which is the typical situation when boundary conditions are not periodic. However, we have 
shown that, for periodic boundary conditions, constant magnetization profiles are locally stable (see Appendix D). 
but the proof of non existence of generic magnetization profiles with larger entropy, and hence the global stability 
analysis, remains to be done. 

An advantage of the method we have exposed is its flexibility and applicability to more complex models. For 
instance, some results of the kind presented here have been already obtained for the a-Blume-Emery-Griffiths model [Sl 
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FIG. 6: Entropy density versus energy density for the a-Ising model but with periodic boundary conditions and with a = 0.8. 
The solid curve represents the theoretical solution, whereas the dashed and dotted lines are determined by microcanonical 
Monte-Carlo simulations with A*' = 34 and A'^ = 100, respectively (data provided by R. Salazar). 

It is also interesting to check how fast one reaches the N —^ oo solution as one increases the number of spins N. 
Fig.^lshows entropy density as a function of energy density for the a-Ising model with a = 0.8 and periodic boundary 
conditions. The figure emphasizes that the asymptotic result is already accurate enough for N = 100. 



CONCLUSIONS AND PERSPECTIVES 



In this paper we have discussed examples of the application of large deviation techniques to the study of the statis- 
tical mechanics of infinite-range models, at equilibrium, in the microcanonical and canonical ensembles. Besides that, 
we have shown how to construct a mean-field Hamiltonian for the Ising model in one dimension with 1/r" interaction 
(0 < a < 1). The solution of simple toy models already shows interesting ensemble inequivalence features, like neg- 
ative specific heat. Among these simple models, one should point out those with continuous state variables, because 
the Hamiltonian dynamics becomes accessible and one can then study also non-equilibrium features. Remarkable is 
also the solution of the Colson-Bonifacio model, which is believed to capture the phenomenology of the saturated 
state. 
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It is important to emphasize that the method proposed in this paper does not apply to all long-range interacting 
systems. In particular, those for which statistical mechanics cannot be reduced to a mean-field variational problem 
are excluded. As presented in Section ^ the method is strongly dependent on the possibility to introduce global 
or coarse-grained variables: examples arc the averaged magnetization, the total kinetic energy, etc. The coarse- 
grained variables allow one to describe structures whose size is of the order of the total size of the system; however 
they may be insufficient to characterize the effect of short-range interaction. A typical example is the Ising model 
with attractive short-range interactions and repulsive long-range couplings studied by Grousson et al. |49j. The 
creation of large scale structures needs an infinite energy, because of the repulsive long-range component. However, 
the competition between repulsion and attraction can lead to interesting structures, such as alternating positive and 
negative magnetized stripes. Indeed, despite the long-range character of the interaction, this system is additive and 
these structures are therefore compatible with the predicted zero magnetization at large scale. Still in the context of 
canonical ensemble, worth quoting are the results of Kardar [53 who, analyzing an Ising model with both short-range 
and long-range interactions was able to derive the exact free energy by a minimization procedure. Obtaining similar 
results for the microcanonical entropy would be extremely important. 

The application of the techniques described in this paper to more realistic A'^-body systems is an important issue. 
The extension to wave-particle interactions [6j should not be too difficult. The mean- field description of the two- 
dimensional point vortices model |5ll | has been rigorously obtained in a series of papers [13. I52L |^ . On the contrary, 
due to the strong short-distance singularity, similar results for self-gravitating systems [2, l54| have not been obtained 
and are a challenging current issue. What is usually done is: i) in equilibrium, to conjecture the validity of the mean- 
field description; ii) out-of-equilibrium, to consider the Vlasov-Poisson equation as a good approximation of the short 
time dynamics. Once the mean- field description is introduced at finite time (Euler equation, Vlasov-Poisson equation), 
the statistical mechanics can be derived using large deviation techniques [iall3|- However, it is not equivalent to the 
original statistical mechanics of the TV-body system. Examples of this inequivalence, due to the exchange of the two 
limits t —> 00 and N —> 00, are given in Refs. [3,|3a| for the HMF and the Colson-Bonifacio model. 
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numerical data used in Fig. Eland A. Alastuey, A. Antoniazzi, P. Cipriani and H. Touchette for discussions. This work 
has been partially supported by the French Ministere de la Recherche grant ACI jeune chercheur-2001 N° 21-31 and 
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Appendix A: Alternative derivation of the microcanonical solution of the HMF model 

The microcanonical solution of the HMF model H32|) can be alternatively obtained using the traditional method 
applied for self-gravitating systems |37l | . Denoting by K and V kinetic and potential energy, respectively, the number 
of microscopic configurations corresponding to the energy E is given by 

^NiE) = fl[dp^de,SiE-HN) (65) 

i 

= Jl[dpideiJdKdlK-Y,^\s{E-K-V{{e,})) (66) 



JdKjl[dp,slK-Y,Y) JUde^HE 

i \ i / i 



-K-V{{e,})) . (67) 



f2ki„(A') 0„,„f(B-A') 

The factor fikin is classical and corresponds to the volume of the hyperspherc with radius R = \/2A' in N dimensions: 
its expression is fikin = 27r^/^/r(l + N/2). Using the asymptotic expression of the F-function, In r(A) ~ A^ In A — A^, 
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one obtains 



i2ki„ (A ) ^ exp ( — 



N 



N 
l + lnTT-ln— +ln(2/\) 



exp ( — [1 + ln(27r) + Inu] 



(68) 
(69) 



where u = 2K/N. Defining the configurational entropy per particle Sconf(V^) = {\'0-Vtconf{NV))/N , where V = 
{E — K)/N , Eq. l|^ can be rewritten as 



ni^{eN) ^ °° — I ducxp 



iV(i + i^ + iln. + w(^) 



Hence, solving the integral in the saddle point approximation, gives the entropy modulo a trivial constant 

S(e) = lim —InnMieN) 



^ u 



lnU + Sconf(V^) 



(70) 

(71) 
(72) 



to which Eq. H43() reduces once the sup on M is performed. Indeed, this expression of the entropy assumes the 
knowledge of the configurational entropy Sconf ■ which is determined by solving the extremal condition in Eq. H42|l . 
and is restricted to total vanishing momentum w = 0. Hence, the method given in the text is slightly more general. 



Appendix B: Mapping the FEL model onto HMF 



The microcanonical solution of the FEL model (See Eq. 1)48(1 ) can be expressed in terms of the HMF Hamiltonian 
using the Laplace representation of the Dirac (5-function and a Gaussian integration. After performing the change of 
variables 9i = 9i — (p, the microcanonical volume of the FEL is given by 



n{E) 



Y[dp,d0,dAS{E- Hn) 



Y\_dpid9idA 



J- fdxAE-HN) 



2i7r 



(73) 
(74) 



where F is a path on the complex A-plane, going from —ioo to +ioo, which crosses the real axis at a positive value. 
Introducing the FEL Hamiltonian 



(75) 




niE)^Jlldp,e \ '-'^JJU d\dA ^JdX e^^-^A^ - ^^^^y) 



and performing the Gaussian integral over the field variable A, one gets 



N 2 
P, 



n{E) 



l[dp, J l[d9.^ J^dX 



.Ml_l 
S 2 



A I S - ^ ^ - iV^ - ^ In [AAr(-<5)] 



N 2 

lldp. I l[d9,S U-^^ -TVll^ - i,\n[XN{-S)] 



3 = 1 



5 2 



(76) 

(77) 
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In the large A^-limit, one can neglect all constants and the InA^ term in the argument of the Dirac S, obtaining the 
microcanonical volume for the Hamiltonian 






(78) 



Hence, for negative values of the parameter 6, solving the microcanonical problem for the FEL Hamiltonian (|48|l is 
formally equivalent to obtaining the solution of the HMF Hamiltonian H33|) with C = 2/6 and M^ = 0. 



Appendix C: Derivation of formula J58I) for the a-Ising model 



Denoting by Ejv ~ {^li 1} the phase space of the a-Ising model, its Hamiltonian is defined as a function of E 



N 



to: 



H 



N 



UJN 



^EV 



bibj 



1->3 



-J\ 



(79) 
(80) 



where lom = (Si, . . . , Sn) is a given microscopic configuration. 

The coarse-graining operator Yj^^k divides the lattice into K boxes of size n ~ N/K and defines a locally averaged 
magnetization rrik in the k box Bk- The coarse-grained magnetization is then a step function, which takes a constant 
value 



rrik 



n ^ — ^ 



iSSfe 



in B/j. It is convenient to introduce the operator 

Yn.k ■■ ^N^L\[Q,l]) 



(81) 



(82) 



which maps the configuration space to the coarse-grained magnetization. The length of the lattice is then renormalized 
to the interval [0, 1] and, in the limit N — > oo, this operator defines a continuous magnetization profile m{x). The 
limit of the number of boxes _ftr -^- oo is taken in such a way that the number of sites per box diverges N/K ^ oo, as 
already mentioned in the text. In the following, we will use free boundary conditions (The choice of periodic boundary 
conditions would change only some details of the calculation) . 

The energy as a function of the coarse-grained magnetization is defined as 



J, 



h[YN,K{i^N)]^ ^ -{l~mkmi)dki , 



where 



d. 



kl 



k,l=l 



k/K pl/K -, 

/ dx dy- 

{k~i)/KJ{i-i)/K F-yr 



(83) 



(84) 



Hence, the composition of the operator Yn^k with the functional h allows to define a map, which associates to each 
microscopic configuration lun G ^n a given energy. Our aim is to find a series K{N) such that 



lim sup 



Hn{ujn) 



7V2-" 



h[YN.KiN){(^N)] 







(85) 
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This is what wc mean by uniform convergence of the Hamihonian Hn{llin) to its functional form H[m(x)], defined 
in formula (|59|l . 

The proof relies on the long-range character of the interaction (0 < a < 1). It is straightforward but lengthy. 

Let us first express h[Yi^ j^ /pf\{u; n)] directly as a function of the spin variables Si 



K 



h\Yj 



J , 



N,K(N) 



MHE3E4(1-^^^: 



n- ^^ 2 



(86) 



In order to compare expression (|86|l with Hn/N'^ ", let us first introduce the reduced Hamiltonian 



9N.K 



1 ^ 

'2-Q Z^ Z^ 



J (1 - S,S,) 



^2-Q Z^ Z^ 2 Infc-n/K 



(87) 



which is in fact the original Hamiltonian Hp^/N'^ " where the distance between two sites is approximated by the 
distance among the boxes to which they belong. Using the triangular inequality, one gets 



Hn{ujn) 



]Sj2-a 



- h[YNJ<{N){'^N)] 



< 



Hn{ujn) 



N^-c 



9n,k 



gN,K - h^N,K(N){^N)\ 



= A + B 



(88) 



We will show how both A and B can be bounded from above by quantities which vanish when N —^ oo, under the 
hypothesis that < a < 1. 

1. Upper bound of A 

Let us first rewrite H^ /N'^~°' as a sum of two terms: the first contains contributions from sites which belong to 
different boxes, the second from those which are in the same box. 
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Then 



A < 
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Ai + A2 .(90) 



The bound of A2 is easy to find since J(l — SiSj)/2 < J and \i — j\" > 1 for a e [0, 1]. One obtains 

A2 < J^ — = J . 

- ^2-a X 



(91) 



This quantity vanishes in the large N limit if K diverges faster than N" (and certainly slower than N). This is the 
first point where the long-range nature of the interaction is used. 

In order to bound Ai let us divide the first sum in formula (|90() into three parts : k — I > 1, k — I < —1 and 
\k — l\ = 1. The last part, \k — l\ = 1, can be bounded from above by 2JN°'/K, which vanishes in the large N limit. 
The two remaining parts are symmetric, hence we treat only the case k — I > 1. Since i G Bk and j G Bi, this implies 
that n{k — / — 1) < i — j < n{k — / -|- 1). Hence, 



1 



< 



< 



n"(fc-/ + l)" \i-j\" n"(fc-;-l)" 



(92) 



Substracting l/\nk — nil", o^^c gets 

1/1 1 



< 



1 1 

< 



n"V(fc-^ + l)" {k-l)"/ \i~j\" n°'{k-l)°' n" \{k - I - 1)°' {k - 1)°" 



This leads to 



Thus 



\i-j\a n'^{k-l)°' 



< 



n" \{k-l- 1)" (fc-/)° 



7V2 



1 ^ 

2-a 2^ 2^ 



k,l=l ieBk 

k-i>ijeBi 



\i-j\a n°'{k-l)°' 



1 n^ ^ 

-a j-,a Z ^ 



2-a K k~2 r 



1 



k-l>l 



—y 

2-a Z-^ 



k=3 1=1 
K 



1 



1 



(fc-/-i)" [k-iy 



K 



fc=3 



1 - 



1 



(fc-l)« 



< 



K 



^2-a 
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(93) 
(94) 

(95) 

(96) 
(97) 



Hence, this part of ^i is bounded by JK"" ^ , which vanishes in the large N limit. The symmetric part of Ai k — l < —1 
can be treated exactly in the same way. 



2. Upper bound of B 



Dividing the sum in two parts, one for sites in different boxes and the other for sites in the same box. and using 
the triangular inequality, one gets 



B < 



N^-c 



N 



J , 



EE^i^-SA 



k,l=l JG-Bfc 
l^k j&Bi 



1 dkl 

''^' y\nk-nl\°' n2 



E§ E |(i-^^^. 



n- ^^ 2 

k ijeSfc 



= Bi+B2 



(98) 



In order to estimate the size of ^2, one has first to evaluate ^^, dkk- Indeed, 



^kl 



1/K fl/K 



JO 



dx dy- 



\x-y + 



k-l I 
K I 



K^ 



du dv- 



\u-v + k-l\°' 



(99) 



Therefore ^^ d^k ~ K°'~^. Hence, B2 < JK°'~^, which vanishes in the large A^-limit. 

For what Bi is concerned, exchanging the modulus with the sums and using the expression for du (formula H99() ). 
one obtains 



TV 



Bi < jJ2 



k,l=l 
l^k 



N'^-°'n°'\k - 11°" 



dki 



< 



w 1 

2^ K'^-a 



k.l=l 

l^k 



/ dudv-. ; 71 7- -T— 



(100) 



Similarly as for Ai, one divides the sum in three parts fc — ^>l,fc — /<— 1 and |fc — Z| = 1. The last part gives a 
term of order JK""^^, which vanishes in the large TV limit. If A; — / > 1, analogously to Eq. H93() . one obtains 



1 i__ r r 1 1^ 1 i__ 

{k-i + 1)« (fc - 0" ^ Jo Jo ""lu-v + k- /|" |fc - ;|" '^ (fc - / - 1)« (fc - ly 



(101) 



which leads to 



du dv 



1 



"'0 



\u-v + k-l\'^ \k~l\^ 



< 



1 



(fc-/-i)" [k-iy 
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(102) 



Summing this term over k and I gives again a factor of order K. Then, this part is also bounded from above by a 
factor JK"^^. The symmetric part A: — ^ < 1 is equally treated. 

All terms in B are bounded by JK°'~^, and therefore vanish in the large TV limit. This concludes the proof of 
formula H85|) and assures that h converges to H. 

We would like to emphasize that throughout of this derivation, the only place where we use the Ising interaction 
of the model is where wc bound J(l — SiSj)/2 by J. It is then easy to realize that the method can be generalized 
to other models like Potts, Blume-Emery-Griffiths, XY, with spatially decaying interactions. On the contrary, the 
presence of the lattice is crucial, because it avoids microscopic configurations where the state variables concentrate 
on single point, leading to a divergence. 

Appendix D: Proof of stability of constant magnetization profiles for the q- Ising model with periodic 

boundary conditions 

We prove here the stability of constant magnetization profiles for the a-Ising model in the context of the canonical 
ensemble. This also implies stability in the microcanonical ensemble. Let us consider the free energy F[m{x)] as a 
functional of m{x) 



F[m{x)] = -S[mix)] + f3H[m{x)] 



(103) 



where the entropy functional iS'[77i(a::)] is defined in formula Hfcil|l . We will limit ourselves to periodic boundary condi- 
tions. Let us study the second variation of F 



S^Fim{x)) 



dx s" [m{x)] STn{x) — P— I dx I dy 



Sm(x)Sm{y) 



(104) 



where s" is the second derivative of s. Let us express the variation 6m{x) in Fourier components 

+ 00 



Sm{x) = Y^ (Jmfce^"'^^ 



fc— — OO 



and define 



(105) 



c^ = / dx 



Jli-nkj. 



k-y|" 



= 2 



^/^ , cos(27rfca;) 
dx r-; 



(106) 



where, since we consider periodic boundary conditions, | • | denotes the distance on the circle. Using the inequality 
s"{ni{x)) < s"(0) and replacing the Fourier expansion (|105|) in Eq. (|104|) . we obtain 



5'F{m.{x))> J2 i-s"{0)-(3-c: 



fe=-c 



J 



(107) 



Remarking that all coefficients c^ are positive reals and that for all fc, c^ < Cg = 2"/(l — a), one can bound from 
above all the coefficients in the Fourier sum H107|l by (— s"(0) — /?Jcg/2). 

Thus for /3 < — 2s"(0)/(Jco) = /3c, one has d'^F{m{x)) > 0. In such a case, the free energy is strictly convex, and 
hence its minimum is unique. This proves the global stability of the state 7n(x) = 0. 
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For f3 > f3c, we do not study global stability. Let us however show the local stability of the uniform magnetization 
state m(x) = m. Using Eq. ()104|l . one has 

6'F{m{x))> ^ l^s"{m)-P-c'^]Sml . (108) 

k=-oo ^ ^ 

Using again the property c^ < Cg, we note that all the Fourier coefficients of Eq. (|1U8|I are bounded from above 
by —s"{m) — /3Jcq /2. Observing that this latter expression is the second variation of the free energy for uniform 
magnetization profiles, i.e. (1 — m^(/3))~^ — (3/ (3c, the analysis of Section HVI implies that this quantity is non negative 
both for m y^ and (3 > [3c and for ttt, = and (3 < Pc- This proves that close to a constant magnetization profile, 
there is no non uniform magnetization profile which gives a smaller free energy functional i^[7Ti(x)]. 
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